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ABSTRACT 

q t We report the discovery of X-ray emission from SN 1970G in M101, 35 yr after its outburst, using deep X-ray 

imaging with the Chandra X-Ray Observatory. The Chandra ACIS spectrum shows that the emission is soft 
{■&! keV) and characteristic of the reverse-shock region. The X-ray luminosity, L 03 _ 2 = (1.1 ± 0.2) x 10 37 ergs 
s' 1 , is likely caused by the interaction of the supernova shock with dense circumstellar matter. If the material was 
deposited by the stellar wind from the progenitor, a mass-loss rate of M = (2.6 ± 0.4) x 10 -5 M e yr' 1 (vj 10 
km s' 1 ) is inferred. Utilizing the high-resolution Chandra ACIS data of SN 1970G and its environment, we 
reconstruct the X-ray lightcurve from previous ROSAT HRI, PSPC, and XMM-Newton EPIC observations, and find 
a best-fit linear rate of decline of L oc t~‘ with index s = 2.7 ± 0.9 over a period of ~20-35 yr after the outburst. 

As the oldest supernova detected in X-rays, SN 1970G allows, for the first time, direct observation of the transition 
from a supernova to its supernova remnant phase. 

Subject headings: circumstellar matter — supernova remnants — supemovae: individual (SN 1970G) — 

X-rays: general — X-rays: individual (SN 1970G, M101) 


Since the launch of the Chandra and XMM-Newton obser- 
vatories, the number of supemovae (SNe) detected in X-rays 
in their near aftermath has more than doubled (see Immler & 
Lewin [2003] for a review). 4 The high-quality X-ray spectra 
have confirmed the validity of the circumstellar interaction 
models (see Fransson et al. 1996, and references therein), which 
predict a hard spectral component (S 10 keV) for the forward- 
shock emission during the early epoch (S 100 days) and a soft 
thermal component (S 1 keV) for the reverse-shock emission 
after the expanding shell has become optically thin. The soft 
emission component dominates the X-ray output of the inter- 
action regions due to its higher emission measure and higher 
electron densities. This expected “softening” of the X-ray spec- 
trum has been observationally confirmed for a number of young 
SNe, such as SNe 1978K (Schlegel et al. 2004), 1979C (Immler 
et al. 2005), 1993J (Zimmermann & Aschenbach 2003), 
1999em, and 1998S (Pooley et al. 2002). 

Where sufficient data are available, X-ray light curves could 
be established over significant timescales (e.g., ~25 yr for SNe 
1978K and 1979C; Schlegel et al. 2004, Immler et al. 2005; 
«8 yr for 1993J; Immler et al. 2001, Zimmermann & Aschenbach 
2003). However, little is known about the transition from a young 
SN into its supernova remnant (SNR) phase, as there have been 
only a few detections of intermediate-age (25-300 yr) SNe in 
the radio and none in X-rays. While the X-ray emission of SNe 
is dominated by the interaction of the shock with the ambient 
circumstellar matter (CSM), likely deposited by the progenitor’s 
stellar wind, the emission from a SNR is thought to originate in 
the shocked interstellar medium (ISM). At an age of «35 yr, SN 
1970G is the oldest SN detected in X-rays and closes this gap, 
allowing one, for the first time, to witness the transition from a 
SN to its SNR phase. 
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As the first SN detected in the radio band (Gottesman et al. 
1972), just 1 month after its peak optical brightness (1970 July 
30; Detre 1970), SN 1970G held the key to our understanding 
of SN events and the interaction of a SN shock with its en- 
vironment. Twenty years after the explosion, SN 1970G was 
rediscovered at 3.5 cm and 20 cm at flux densities similar to 
other radio SNe, such as SNe 1950B, 1957D, 1961V and 1968D 
at the same epochs in their evolution, but showed a remarkably 
fast rate of decline, with a power-law index of —1.95 ± 0.17 
(Cowan et al. 1991). Subsequent radio monitoring confirmed 
this fast rate of decline, followed by a flattening of the radio 
light curve at an epoch of «30 yr, which was discussed in the 
context of the SN shock running into the denser ISM, indicating 
the onset of the SNR phase (Stockdale et al. 2001). 

In this paper we report on the Chandra X-ray observation 
of SN 1970G, as well as previous X-ray data from ROSAT and 
XMM-Newton. In § 2 we briefly describe the data and analysis 
thereof and report on the X-ray spectrum and multimission X- 
ray light curve in § 3. We discuss the results in the context of 
the CSM interaction model and put SN 1970G in the context 
of the evolution from a SN to a SNR in § 4, followed by a 
summary in § 5. 

2. DATA PROCESSING AND ANALYSIS 

SN 1970G has been observed with the Chandra X-Ray Ob- 
servatory as part of the 1 Ms observation of the host galaxy 
M101 (NGC 5457; K. D. Kuntz 2005, in preparation). Five <p 
individual Chandra ACIS observations (sequence numbers 
600389, 600389, 600389, 600389, 600389) from 2004 July 5- 
1 1 were merged into a single observation, resulting in a cleaned 
and exposure-corrected on-source exposure time of 139.5 ks 
for the ACIS-S2 chip. Details of the data calibration and anal- 
ysis will be explained in K. D. Kuntz (2005, in preparation). q 4 
The high-resolution (=0'.'5 FWHM) Chandra data were used 
to search for X-ray emission from the position of SN 1970G, 
and to examine the contamination by nearby X-ray sources. A 
Chandra ACIS-S2 image of SN 1970G and its environment is 
given in Figure 1. q s 

We further constructed the long-term X-ray light curve of 
SN 1970G by extracting source counts from previous ROSAT 
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Fig. 1 . — Full-band (0.3-8 keV) Chandra ACIS-S2 image of the region (size 
lxl) around the position of SN 1970G. The image is smoothed with a 
Gaussian of 2”5 (FWHM) and is given in logarithmic gray scale ranging from 
0-1.5 counts pixel' 1 (1 pixel corresponds to 0'.'5). Contour lines are 0.2, 0.4 
(black), 0.6, 0.8, 1, 1.2, 1.4, and 1.6 counts pixel" 1 (white). The radio position 
of SN 1970G (Cowan et al. 1991) northeast of the H ti region NGC 5455 is 
marked by a white cross. 

HRI (sequence numbers RH6OO820NOO, RH600820A01; 
merged exposure time 176.1 ks) and ROSAT PSPC observa- 
tions (sequence number RP600108NO0; exposure 34.5 ks). Ar- 
chival XMM-Newton EPIC observations (ObsID 0104260101, 
43.3 ks; 0164560701, 41.8 ks; 0212480201, 32.4 ks) were also 
analyzed according to standard analysis procedures described 
in the XMM-Newton ABC Guide, version 2. 01. 5 

A short (2.8 ks) Einstein HRI observation 154 days after the 
outburst of SN 1970G did not result in the detection of the SN 
and is not included in this study. Archival ASCA observations 
were not used due to the large aperture («1' FWHM) and 
contamination by nearby X-ray sources. The ROSAT a.nd XMM- 
Newton count rates were further compared to the published 
results by Wang et al. (1999; ROSAT) and Jenkins et al. (2005; 
XMM-Newton). 

3. RESULTS 

A pointlike X-ray source is detected at the radio position of 
SN 1970G (Cowan et al. 1991), with an offset (S0"5) similar 
to the spatial resolution and astrometry of the Chandra obser- 
vation (see Fig. 1). SN 1970G is spatially resolved and well 
separated from a nearby (=5") H n region NGC 5455. Extract- 
ing exposure-, background-, and aperture-corrected counts from 
the position of SN 1970G gives a count rate of (4.3 ± 
0.6) x 10'“ counts s' 1 . While the limited photon statistics do 
not allow a detailed characterization of the energy distribution 
of the recorded photons, spectral analysis shows that the emis- 
sion is soft, with the bulk of the emission being confined to 
the -&2 keV band. 

Adopting a kT = 0.6 keV thermal plasma spectrum typical 
for the late-time (felOO days) emission of SNe and assuming 

5 Available from http://heasarc.gsfc nasa gov/docs/xmm/abc/. 



Fig. 2. — Soft-band (0.3-2 keV) X-ray luminosities of all SNe detected to 
date (filled squares) and historical SNRs (filled triangles) as a function of age 
(in units of years). The X-ray light curve of SN 1970G is marked by filled 
diamonds with error bars (left to right'. ROSATPSPC, ROSATHRI and Chandra 
ACIS). XMM-Newton EPIC upper limits are indicated by arrows. 

a Galactic foreground column density with no intrinsic ab- 
sorption ( N h = 1.16 x 10 20 cm' 2 ; Dickey & Lockman 1990) 
gives a 0.3-2 keV flux and luminosity of / 03 _, = (1.8 ± 
0.3) x 10' 15 ergs cm' 2 s' 1 and L 03 _ 2 = (1.1 ± 6.2) x 10 37 
ergs s' 1 , respectively, for a distance of 7.2 Mpc (Stetson et al. 

1998) . 

Since no information is available about the spectral evolution 
of SN 1970G or other SNe at this age, we used the same spectral 
template to convert ROSAT and XMM-Newton count rates into 
fluxes and luminosities to construct the long-term («20-35 yr) 
X-ray light curve. 

An X-ray source is visible in the ROSAT HRI images of 
M101 at the position of SN 1970G (see Fig. 2 in Wang et al. 

1999) . The source, however, is not included in the ROSAT HRI 
catalog since it is slightly below the employed detection thresh- 
old (S/N = 3.5; Wang et al. 1999). Extraction of source counts 
within the HRI aperture and subtraction of emission from the 
nearby H ii region NGC 5455 (L 0 3 _ 2 = 2.2 x 10 36 ergs s' 1 ) 
gives a luminosity of L 03 _ 2 = (2.6 ± 0.6) x 10 37 ergs s' 1 . 

An X-ray source is also listed in the ROSAT PSPC catalog 
at the position of SN 1970G (source P12; Wang et al. 1999), 
with a signal-to-noise ratio of S/N = 4.5. We used the high- 
resolution Chandra observation and subtracted the integrated 
flux of all detected Chandra sources from the ROSAT PSPC 
flux within the PSPC aperture («25" FWHM). Assuming that 
the residual emission arises from the SN itself, we obtain a 
luminosity of L 03 _ 2 = (4.9 ± 1.1) x 10 37 ergs s' 1 . 

A faint X-ray source is visible in each of the nine XMM- 
Newton pn, MOS 1, and MOS2 images at low energies (<2 keV). 
No photons were recorded in the hard (2-6 keV) band. The 
source, however, is extended («20") and shows no concentration 
of the recorded photons at the radio position of SN 1970G. Since 
it is likely that the emission is contaminated by the nearby H n 
region and other soft emission components inside the Galactic 
disk due to the larger aperture of the EPIC instruments (~6" 
FWHM), we only use the inferred count rates as upper limits. 
Discrepancies between the published pn and MOS count rates 
for source 44 (factor of 16.5; Jenkins et al. 2005), corresponding 
to SN 1970G, for one of the published XMM-Newton data 
(ObsID 0104260101) make a comparison problematic. 

We calculated the mass-loss rate of the progenitor as a func- 
tion of the stellar wind age using the relationship L x = 
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TABLE 1 

X-Ray Properties of SN 1970G 


Epoch Count Rate f 0 j_, L J 3 _ 2 M 

(yr) Date Instrument (1 O' 4 counts s' 1 ) (10" 15 ergs cm -2 s' 1 ) (10“ 37 ergs s" 1 ) (10 _i M e yr' 1 ) 


20.9 1991 Jun 8-9 ROSAT PSPC 9.0 ± 2.1 7.9 ± 1.8 4.9 ± 1.1 3.2 ± 0.7 

26.1 1996 May 14-Nov 23 ROSAT HRI 1.8 ± 0.4 4.2 ± 0.9 2.6 ± 0.6 2.7 ± 0.6 

31.9 2002 Jun 4 XMM-Newton EPIC <48.3 <10.2 <6.3 <5.8 

34.0 2004 Jul 5-11 Chandra ACIS 4.3 ±0.6 1.8 ± 0.3 1.1 ± 0.2 2.6 ± 0.4 

34.0 2004 Jul 23 XMM-Newton EPIC <80.6 <17.0 <10.6 <7.7 

34.5 2005 Jan 8 XMM-Newton EPIC <119.3 <25.2 <15.6 <9.4 


4/(7r/n 2 )A(T) x ( Mlv w ) 2 x (v s t)~\ where m is the mean mass 
per particle (2.1 x 10“ 27 kg for a H+He plasma), A(T) the 
cooling function of the heated plasma at temperature T, M the 
mass-loss rate of the progenitor, v„ the speed of the stellar wind 
blown off by the progenitor, and v, the speed of the outgoing 
shock (see Immler et al. 2002). A constant shock velocity of 
v s = 9000 km s" 1 , similar to other core-collapse SNe, such as 
SN 1979C (Marcaide et al. 2002) was assumed. An effective 
(0.3-2 keV band) cooling function of A = 3 x 10“ 23 ergs 
cm" 3 s" 1 for an optically thin thermal plasma with a temperature 
of 10 7 K was adopted, which corresponds to the assumed ther- 
mal plasma temperature. Assuming different plasma temper- 
atures in the range 0.5-1 keV would lead to changes in the 
emission measure of £10%. Key observational properties of 
SN 1970G are listed in Table 1. 

4. DISCUSSION 

Previous X-ray observations of SN 1970G with ROSAT , 
ASCA, and XMM-Newton lacked the spatial resolution needed 
to separate the SN from the nearby (=5") H n region NGC 
5455. Using the subarcsecond Chandra imaging capability, in 
combination with the deep exposure time (140 ks), we detect 
soft (£2 keV) X-ray emission from SN 1970G with a lumi- 
nosity of L 03 _ 2 a 1 x 10 37 ergs s" 1 , 35 yr after its outburst. 
The source is spatially resolved from the nearby H n region 
NGC 5455 (=5"), which is observed at a £20% flux level 
compared to SN 1970G. 

Given the spatial information from the high-resolution Chan- 
dra imaging, we recover X-rays from SN 1970G in the ROSAT 
HRI and PSPC observations. SN 1970G is not conclusively 
detected in three recent (2002-2005) XMM-Newton observa- 
tions due to the extent of the source (=20") as imaged by the 
pn and MOS instruments, and positional offsets larger than the 
instrumental uncertainties (=5"). 

In combination, the multimission X-ray data show a constant 
rate of decline which confirms that the emission is due to the 
SN itself. The rate of decline of L oc t~ s with index s = 
2.7 ± 0.9 is rather steep and similar to the X-ray rates of de- 
cline observed for the Type Iln SNe 1986J (L x oc r~ 3 ; Temple 
et al. 2005) and 1988Z (f 39 ; Aretxaga et al. 1999). Since we 
cannot conclusively exclude additional contamination of the 
earlier ROSAT data with unresolved X-ray emission compo- 
nents in M101, our inferred rate of decline is rather 
conservative. 

The rate of decline inferred from the X-ray data is consistent 
with the slope of the radio decay (—1.95 ± 0.17; Stockdale et 
al. 1991), which confirms that the emission probes similar 
regions in the shocked CSM. The X-ray data, however, do not 
show the flattening of the decay at an age of =30 yr, as seen 
in the 20 cm radio data (—0.29 ± 0.13). While this might be 
due to the limited number of X-ray data points, it could also 
indicate that the reverse shock, which produces the X-rays and 


travels at a speed of some 1000 km s" 1 slower than the forward 
shock, has not yet reached the dense ISM, which is probed by 
the radio emission originating in the forward shock. 

The lack of X-rays above 2 keV confirms that the emission 
is soft and likely thermal. Such a soft spectrum is expected for 
the late emission of a SN, where the shocked plasma emission 
is the product of the reverse shock interaction. 

Assuming that the emission arises from shock heated plasma 
deposited by the progenitor’s stellar wind, a mass-loss rate of 
M = (2.6 ± 0.4) x 10" 5 M @ yr" 1 (vJlO km s" 1 ) is inferred 
for the late epoch observed with Chandra. The mass-loss rates 
inferred from the ROSAT PSPC and HRI data [(3.2 ± 0.7) and 
(2.7 ± 0.6) x 10" 5 M q yr" 1 (v w /lQ km s" 1 )] are consistent 
with a constant mass-loss rate over a period of =20-35 yr after 
the peak optical brightness. This corresponds to =19,000- 
31,000 yr in the stellar wind age, using the conversion t„ = 
tvjv n and assuming a stellar wind speed of 10 km s" 1 and a 
shock front velocity of 9000 km s" 1 . 

Our inferred mass-loss rate is consistent with that derived 
from radio monitoring [2 x 10“ 5 M 0 yr" 1 (ij/10 km s" 1 ); 
Weiler 1993]. Furthermore, the measured mass-loss rate for SN 
1970G is similar to those inferred for other young (£20 yr) 
Type II SNe, which typically range from 10" 5 to 10' 4 M q yr" 1 
(see Immler & Lewin 2003). This indicates that the X-ray 
emission arises from shock-heated CSM deposited by the pro- 
genitor rather than shock-heated ISM, even at this late epoch 
after the outburst. 

In the absence of a dense ISM, the X-ray evolution of a 
core-collapse SN into a SNR can be described by a simple 
power law. X-ray luminous SNe, such as SN 1993J, can there- 
fore evolve into X-ray-luminous SNRs such as Cas A 
(L 0} _ 2 = 2.8 x 10 37 ergs s" 1 at an age of =320 yr) with a 
swept-up gas mass around 4 M 0 (Dunne et al. 2003) if the 
measured mass-loss rate of =10" 4 M @ yr" 1 is sustained over a 
period of some 10 4 yr. 

In Figure 2 we plot the X-ray evolution of all SNe detected 
to date, 6 as well as the X-ray luminosities of historical SNRs in 
the Galaxy, LMC, SMC (as observed with Chandra 7 ), M31 (Sup- 
per et al. 2001; Kong et al. 2003), and M33 (Ghavamian et al. 
2005). A comparison of the X-ray luminosities of the core-col- 
lapse SNe (visible in the left part of Fig. 2) with well-evolved 
SNRs (right part of Fig. 2) shows that the transition from a SN 
to a SNR appears to be rather smooth (see also Fig. 2 in Stockdale 
et al. 2001). If the evolution of SN 1970G continues as observed 
to date, it will be similar to X-ray-faint Galactic SNRs such as 
SNRs 1181 and 1006 (L x = 10 34 ergs s" 1 ). If the light curve 
flattens significantly after encountering a dense ISM (as in the 
case of SNRs as a result from Type la SNe such as Tycho), SN 
1970G might evolve into a more luminous SNR such as Kepler 
or the Crab (=10 36 ergs s" 1 ). The observed flattening of the radio 

s See http://lheawww.gsfc.nasa.gov/users/immler/supernovaeJist.html. 

7 See http://www.astro.psu.edu/users/green/Main/main5.html. 
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q7 light curve with a power-law slope change from -1.95 ± 

0.17 to -0.29 ± 0.13 at an age of =30 yr (Stockdale et al. 1991) 
might indicate such as scenario. 

X-ray-bright SNe such as SNe 1978K (Schlegel et al. 2004) 
and 1979C (Immler et al. 2005), which show no signs of an 
X-ray decline out to large distances (S10 17 cm) from the sites 
of the explosion, can evolve into luminous SNRs such as Cas 
A through their strong interactions with dense stellar winds 
alone. In case the shock encounters a dense ISM at a later 
stage, these SNe have the potential to attain even higher X-ray 
luminosities as observed for the brightest SNRs (such as W49B 
qs or N132D). 

Clearly, long-term monitoring of known X-ray SNe and new 
detections of intermediate-age (around 100 yr) SNe in X-rays 
are needed to address these questions and to study the transition 
from a SN into a SNR in more detail. Deep galaxy surveys 
with Chandra and XMM-Newton and future X-ray missions 
such as Constellation-X and XEUS will provide the capabilities 
to detect these intermediate-age SNe in the near future and to 
fill the void in the evolution of a SN to a SNR. 

5. SUMMARY 

Key results based on our deep (140 ks) X-ray observation 
of SN 1970G in M101 are: 

1. SN 1970G is detected with a luminosity of L 0J _ 2 = 
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(1.1 ± 0.2) x 10 38 ergs s~\ 35 yr after its outburst in our deep 
Chandra observation of M 101. 

2. SN 1970G is recovered in previous ROSAT HRI and 
PSPC observations and shows a best-fit linear rate of decline 
of L oc t~* with index s = 2.7 ± 0.9 over the observed period 
of 20-35 yr after its outburst, consistent with the observed 
radio decay for an epoch S30 yr. 

3. The emission of SN 1970G is soft (£2 keV) and indicates 
that it originates in the reverse-shocked region, as expected for 
the late emission in the evolution of a SN. 

4. A mass-loss rate of M = (2.6 ± 0.4) x 10~ 5 M 0 yr -1 
is inferred consistent with being constant over a period of 
19,000-31,000 yr in the evolution of the progenitor. 

5. At an age of 35 yr, the SN shock has not yet encountered 
a denser ISM, which would mark the start of the SNR phase 
and is likely to occur at a later stage in the evolution (=50- 
100 yr). 

6. Comparison of the X-ray light curves of all SNe detected 
to date with the X-ray luminosities of SNRs indicates a smooth 
transition of SNe into their SNR phases. 
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